Invalidation of the Kelvin Force in Ferrofluids 
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Direct and unambiguous experimental evidence for the magnetic force density being of the form 
M^VB in a certain geometry — rather than being the Kelvin force MVif — is provided for the 
first time. (M is the magnetization, H the field, and B the flux density.) 



In studying polarizable and magnetizable condensed 
systems, one of a few key concepts is the force exerted by 
applied electric and magnetic fields. A clear and precise 
understanding here is therefore of crucial importance. 
For a neutral, dielectric body, this force is usually taken 
to be the Kelvin force [pi 0, HI, 



(1) 



where P is the polarization, M the magnetization, E 
the electric and H the magnetic field; //q and eo are as 
usual the vacuum permeability and permittivity, and the 
integral is over the volume of the body. (Summation over 
repeated indices is always implied.) 

In spite of general believes, however, the theoretical 
foundation for this expression is far from rock solid, and 
recent considerations have pointed to the validity, under 
conditions to be specified below, of the variant expres- 
sion 01 



[(P,/e )VA + M,V^]d 3 r, 



(2) 



where D e E + P and B /i (H + M) are the di- 
electric displacement and magnetic flux density, respec- 
tively. This circumstance makes a direct measurement of 
the electromagnetic force desirable. 

In this Letter, such an experiment for the magnetic 
part of the force is reported. It consists of a pendulum 
exposed to a well defined, horizontal magnetic field with 
a gradient parallel to the field, see Fig. 1. The weight 
of the pendulum is a disk-shaped container. It is filled 
with ferrofluid and has its face aligned normal to the 
field. The disk is drawn by the field gradient, but leaves 
the applied field essentially unchanged. Given the careful 
design of the experiment, with especially the string of the 
pendulum close to 8 m, its displacement yields an easy 
and accurate measurement of the magnetic force. For the 
given geometry, Eq (||) was verified to great precision. 

The quest for the correct form of the electromagnetic 
force sounds, and is indeed, vintage 19 th century. So 
it is truly amazing that something as elementary as the 
rest position of a pendulum subject to a static magnetic 
field cannot be predicted in general consensus even today. 
Kelvin was the first to have made a lasting contribution 
in this context by first calculating, microscopically, the 
force exerted by an external field on a single dipole, and 
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FIG. 1: Schematic representation of the experiment, a pen- 
dulum deflected by a magnetic force. The inset shows the 
two magnetic coils producing the field. The object marked as 
pendulum is a disk containing ferrofluid, with its face perpen- 
dicular to both the field and field gradient. 



then taking the total force as the sum of the forces on 
all dipoles. Clearly, this is only justified if the dipoles 
are too far apart to interact, if the system is electromag- 
netically dilute. This implies a small magnetization and 
susceptibility, M <C H, x ^ 1- (We shall only be dis- 
cussing the magnetic force from here on.) Part of the 
terms neglected must be quadratic in M, as they nat- 
urally account for interaction. Therefore, this approach 
cannot distinguish between Eq (|]) and (||) , the difference 
of which is |/ioVM 2 , a term quadratic in M, or x- 

The classic derivation of the electromagnetic force is 
thermodynamic in nature, and rather more generally 
valid, cf §15 and 35 of the book on macroscopic elec- 
trodynamics by Landau and Lifshitz 0. Starting from 
the Maxwell stress tensor n^, this derivation subtracts 
the hydrostatic "zero-field pressure" - the pressure that 
would be there in the absence of fields, for given tem- 
perature and density - and identifies the gradient of the 
rest as the electromagnetic (or Helmholz) force. If the 
susceptibility x ~ irrespective of its value - is taken to 
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FIG. 2: Magnetization of a typical ferrofluid (No 7 in Table 2) 
as a function of the field H. The inset shows the dependence 
of x on the concentration of magnetic particles for four fluids 
of similar microscopic make-up, especially the mean particle 
diameter. 



be proportional to the density (or concentration of mag- 
netic particles in a suspension such as ferrofluid), this 
expression reduces to the Kelvin force, Eq (Q). With the 
difference between Eq (|l|) and (||) appreciable if % > 1, 
this derivation may be taken to clearly rule out Eq (||). 

Ferrofluids are stable suspensions of magnetic nano- 
particles, with an initial susceptibilities up to about 5 
that are usually proportional to the volume concentration 
pi of magnetic particles at small fields, see Fig 2. It 
is therefore standard assumption that the Kelvin force, 
Eq (|l|), may be employed wherever x ^ Pi holds ||. 

Closer examination 0, however, has raised doubts 
about the validity of Eq (|l]) for highly magnetizable sys- 
tems, by unearthing an implicit assumption of small sus- 
ceptibilities, made during the derivation as given by Lan- 
dau and Lifshitz. So, in spite of appearance, neither 
does this derivation distinguish between Eq (|l|) and (|^). 
Rather more important, however, is the fact that this 
force is not operative in equilibrium, when the divergence 
of the Maxwell stress vanishes, VjHj = 0, and the gra- 
dient of the zero field pressure compensates the Kelvin 
force. 

Generally speaking, while the electromagnetic force 
is a simple and unique quantity in microscopic physics, 
given by the Lorentz force, it is a multi-faceted, heuristic 
concept in macroscopic physics. Not aware of this fact, 
and possibly for want of something better, we tend to 
take Eq (Q) - with its curious preference of E, H over 
D, B - as the electromagnetic force. (There is of course 
also the macroscopic Lorentz force, which is however zero 
if the system is neutral and insulating). 

To understand why this is a conceptual pitfall, consider 
the hydrodynamic theory for ordinary fluids, such as wa- 
ter or air, in the absence of fields. It consists of a set of 
partial differential equations including the Navier- Stokes 
equation |?]] . By solving these equations with the appro- 



priate initial and boundary conditions, one is in princi- 
ple able to predict any experimental outcome - say the 
trajectory of an airplane - without ever the need to in- 
troduce the concept of force. For the convenience of rea- 
soning and arguing, however, we do label certain terms 
as force densities: in the differential equations, boundary 
conditions, or the solution under consideration. As a re- 
sult, these forces depend on the context and geometry - 
take the form of the airfoil and the associated lift. 

Circumstances are similar in the hydrodynamic theory 
of ponderable systems |l|, ^, and there is just as little 
reason to expect the existence of a unique electromag- 
netic force that is independent from geometry and con- 
text. But we do also need the concept of force: Physics 
does not comprise solely of mathematics, and clear think- 
ing about the valid expressions for the electromagnetic 
force under given conditions will give us considerable 
heuristic power in prediction — without having to solve 
the set of partial differential equations each and every 
time. A simple question about the rest position of a pen- 
dulum subject to a field should be answerable by consid- 
ering an elementary force equilibrium, between the strain 
in the string, the gravitational and the electromagnetic 
force. Solving the hydrodynamic theory with appropriate 
boundary conditions ji) and assuming linear constitutive 
relation, the magnetic force is found to be 



X £n7(l + x)HdA, (3) 



where the surface integral is to be taken over the sur- 
face of the magnetizable body. The subscripts t and n 
denote the tangential and normal field component with 
respect to the surface. As the above fields are continu- 
ous across the interface, they may be both the external 
or the internal field. Eq (||) is rather generally valid. It 
holds irrespective of x' s si ze ? or its dependency on ei- 
ther the density or the volume concentration of magnetic 
particles; it is also easily generalized for nonlinear con- 
stitutive relations. Equilibrium, however, is an essential 
requirement. 

Employing the Gauss law, and provided x is spatially 
constant, with M = xH = x^/mo(1+x)? we mav rewrite 
Eq (||) as a volume integral, 



T = J (/jLoMtVHt + M n VB n ) d 3 r. 



(4) 



If the field is normal to the surface, as in the above ex- 
periment, the first term vanishes, and T w of Eq (^) is 
retrieved; if the field is tangential to the surface, Eq (|l|) 
and Tk hold. 

Casting doubts - however well theoretically founded 
- on the applicability of the Kelvin force certainly re- 
quires experimental validation. Fortunately, in strongly 
magnetizable systems such as ferrofluids, the difference 
between Ty, and T w is appreciable. For linear consti- 
tutive relations, the latter is larger by the factor 1 + \\ 
With the magnetic field perpendicular to the disk sur- 
face, we have B in = B ex = n$H e ^ where the subscripts 
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i n and ex denote the internal and external field, respec- 
tively. This implies H in = if ex /(l + x) an d M = X^in 
= xH GX /(l + x). As the disk containing the ferrofluid is 
not infinite, the external field H ex is not quite the applied 
field Ho. Approximating the disk as a (flat) ellipsoid, we 
may write H ex / (1 + x) = ^in = H$/(l + Dx) , where D is 
the demagnetization factor. (The ferrofluid sample has 
the form of a flat disk with a diameter of a = 50 mm, 
thickness c = 1 mm, and b = a/c = 50, from which 
the demagnetization factor in z direction is calculated as 
D = 0.9694, employing the approximation H D = 1 — 
n/2b +2b 2 .) Inserting these formulas into Eqs (|l|) and 
(0), respectively, we find 

Tk = [^ M oX/(l + Dx) 2 } VHl (5) 
T Y = (1 + x)[^ M ox/(l + D X f]VHl (6) 

They show T w = (1+x)«^*k, and give the force in terms of 
simple, measurable quantities: susceptibility x? applied 
field i^o, the demagnetization factor D, and the ferrofluid 
volume V. (Due to the small gradient of the applied field, 
the volume integral was approximated by multiplication 
with V.) 



description 


symbol 


value 


length of pendulum 


1 


7650mm 


diameter of fluid disk 


a 


50mm 


thickness of fluid disk 


c 


1mm 


volume of fluid disk 


V 


2ml 


mass of pendulum 


m 


35g 


demagnetization factor 


D 


0.9694 


field parameters 


OL\ 


2.5 x 10 4 m- 3 






0.47m 




as 


750m- 1 



Table 1. Technical data of the experiment 



The magnetic field, with its direction and gradient 
both parallel to the ^-direction, is produced by an ar- 
rangement of two coils, which is simply half a Fanselau- 
arrangement JnJ, see Fig. 1. The field is homogeneous in 
the xy-pldine to an accuracy of better than 0.1 %. Both 
the field and its gradient are proportional to the applied 
current /. The field strength shows a quadratic decrease 
in z-direction, leading to a linear dependence of the field 
gradient. Choosing the center of the larger coil as the ori- 
gin of the coordinate system (z = is at the outer edge 
of the coil), the field is parameterized with three coeffi- 
cients, ai — 2.5 x 10 4 m -3 , a 2 = 0.47m, as = 750m _1 , 
as 

H = [ai(z - a 2 ) 2 - a 3 ]I. (7) 

(Ho denotes the z-component of the field. The other com- 
ponents are not zero for z ^ 0, as V • B = 0. The force 
density from the x, ^-components of the field is quadrat- 
ically small, and integrates to zero for a disk centered at 
z = 0.) 



The disk containing the ferrofluid is attached to two 
nylon strings, of length I = 7650 mm and mounted to 
a tripod, forming a pendulum. The rest position of the 
pendulum is at x : y = and z = 7.6 cm. If a mag- 
netic field is applied, the magnetic force acting on the 
fluid attracts the disk towards the coils. This results 
in its displacement and, as shown in Fig. 1, a restoring 
gravitational force, the relevant component of which is 
mg tan(/? = mgAz/l, with m the mass of the disk, g the 
gravitational acceleration, ip the angle of deflection, and 
Az the displacement. In equilibrium, the gravitational 
force is equal to the magnetic one, either Tk of Eq (||) or 
of Eq (^) , but always ~ Vi^o • Writing the magnetic 
force as /3\7Hq and equating it with mgAz/l, we find 

Az = {pi/mg)VHl = XH VH . (8) 

The slope A = 2/31 /mg is measured by plotting Az over 

Due to the great length of the pendulum, even small 
changes of the magnetic force will provide easily mea- 
surable changes of Az. For instance, a change of the 
magnetic force of 10 _5 N corresponds to Az = 1mm. In 
addition, since Az = 13 mm corresponds to cp = 0.1°, 
even relatively large displacements, of the order of a few 
cm, will not lead to a tilt of the disk relative to the field 
direction. This ensures the homogeneity, over the sam- 
ple, of the magnetic field Ho in x-direction. The position 
of the sample is monitored by a digital video system, al- 
lowing a determination of the displacement Az with an 
accuracy of approximately 0.2 mm. 



ferrofluid 


X 


concentration 


particle size 


trade name 


1 


0.183 


1.8 vol% 


7.1 nm 


APGS20n 


2 


0.393 


1.8 vol% 


9.2 nm 


APG510A 


3 


0.543 


1.4 vol% 


9.2 nm 




4 


0.723 


3.7 vol% 


9.2 nm 


APG511A 


5 


1.141 


5.4 vol% 


9.1 nm 


APG512A 


6 


1.27 


6.2 vol% 


8.7 nm 




7 


1.57 


7.3 vol% 


9.2 nm 


APG513A 


8 


1.84 


7.3 vol% 


8.9 nm 




9 


2.27 


7.3 vol% 


10.4 nm 


EMG905 


10 


3.65 


13 vol% 


8.5 nm 




11 


3.95 


16.4 vol% 


9.5 nm 


EMG900 



Table 2. Data of the employed ferrofluids 



With the experimental setup described above, we have 
carried out a series of experiments using 11 different 
ferrofluids having different initial susceptibility, ranging 
from 0.18 to approximately 4, cf table 2. (For ferroflu- 
ids supplied by Ferrofluidics, the trade names are given 
in table 2.) Because only weak magnetic fields up to a 
maximum of Hq = lOkA/m are used in the experiment, 
the validity of linear constitutive relation is ensured, cf 
Fig 2, and we can use the initial susceptibility of table 2 
throughout the evaluation of the data. 
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FIG. 3: Dependence of the displacement Az of the pendulum 
on HqVHq, for fluids 2, 7 & 11, with linear fits. 



Fig. 3 shows the displacement Az of the disk as func- 
tions of HqVHq, for three different fluids from table 2. 
It is obvious that the linear relation between Az and 
Hq\7Ho is well satisfied over the whole range of inves- 
tigation, attesting to the validity of the linear constitu- 
tive relation employed in the evaluation. From plots like 
these, we can determine the slope to an accuracy of about 
3%. By averaging over five runs for each fluid, the accu- 
racy is increased to approximately 1.5%, demonstrating 
the reproducibility of the data. 



Contact with theory is made by plotting A as a function 
of x f° r a ll 11 fluids described in table 2, see Fig. 4. The 
solid line gives A as calculated with T w of Eq (o) , while 
the dashed line is calculated with Tk of Eq (||) - both 
using the data of table 1 and the independently measured 
susceptibility of table 2. 
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Clearly, the agreement between the solid line and the 
experimental data is excellent. Moreover, it is obvious 
that the error committed by using the Kelvin force in- 
discriminately can be considerable. As the correction is 
1 + X? however, this is an issue only for highly magnetiz- 
able systems, say with susceptibilities > 0.2, which helps 
to explain why this was not seen earlier. 



FIG. 4: A = Az/(HoVHo) as a function of the susceptibility. 
The dashed line is calculated employing the Kelvin force Fk, 
of Eq (|H|), while the solid line is calculated with the variant 
form of Eq (|). 
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